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Abstract 

Scaling studies of inclusive quasielastic electron scattering reactions have been used in the past 
as a basic tool to obtain information on the nucleon momentum distribution in nuclei. However, 
the connection between the scaling function, extracted from the analysis of cross section data, 
and the spectral function only exists assuming very restricted approximations. We revisit the basic 
expressions involved in scaling studies and how they can be linked to the nucleon momentum distri- 
bution. In particular, the analysis applied in the past to the so-called scaling region, i.e., negative 
values of the scaling variable y, is extended here to positive y, since a "universal" superscaling 
function has been extracted from the analysis of the separated longitudinal data. This leads to 
results that clearly differ from the ones based solely on the negative-y scaling region, providing 
new information on how the energy and momentum are distributed in the spectral function. 
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I. INTRODUCTION: BASIC ASPECTS OF SCALING 



Scaling studies of inclusive quasielastic (QE) electron-nucleus scattering have largely been 
considered to provide a powerful tool for extracting the momentum distribution of nucleons 
inside nuclei [l|-|7j|. Such analyses have been applied to few-body systems, complex nuclei 
and nuclear matter with an important effort devoted to estimating binding corrections, 
and in particular, the high-momentum components of the nucleon momentum distribution 
which are governed by Short Range Correlations (SRC) j^l, 9]. However, caution should be 
borne in mind for the conclusions reached, since a close relationship between the momentum 
distribution and the scaling function only emerges after some approximations are made. 
These are linked not only to the general description of the electron scattering reaction 
mechanism, but also to the integration limits involved and the behavior of the spectral 

n 

function 

The phenomenon of y-scaling emerges from the analysis of QE (e, e') reactions. The 
scaling function, defined as the QE (e, e') differential cross section divided by an appropriate 
factor involving the single- nucleon cross section jliliol 12 1. is shown to depend only on a single 
variable, y, given as a particular combination of the two independent variables in the process, 
namely the energy and momentum transfers, u and q. In the QE domain and for values 
of u and q large enough, the basic mechanism in (e, e') reactions on nuclei corresponds to 
elastic scattering from individual nucleons in the nuclear medium with "quasi-free" ejection 
of a nucleon from the nuclear system. This implies that the inclusive (e, e') cross section 
is mainly constructed from the exclusive (e, e'N) process, including the contribution of all 
nucleons in the target and integrating over all (unobserved) ejected nucleon variables. In 
other words, QE scattering off a nucleus is simply described as an incoherent sum of single- 
nucleon scattering processes. This approach, which constitutes the basis of the Impulse 
Approximation (IA), although being an over-simplified description of (e, e') reactions, has 
demonstrated its validity under appropriate kinematic conditions. Mechanisms beyond the 
IA (correlations, Meson Exchange Currents (MEC), re-scattering processes, etc.) may play a 
significant role in electron scattering, and hence may lead to non-negligible scaling violations. 

The IA provides an intuitive explanation on how the scaling behavior emerges from the 
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analysis of data. In this case the QE (e, e') cross section is given by 



da 



de'dQ' 



(e,e') 



A 

E 

i=l 



pdpdS 



E 



QP 2 Nl 



da 



de'dQ'dpNidflNi 



(e,e'Ni) 



where the sum extends to all nucleons in the target and {e', Q'} refer to the scattered electron 
variables. The integration over the ejected (unobserved) nucleon variables {p^, E^^fl^} 
has been expressed in terms of the residual nucleus' excitation energy £ and the missing 
momentum p. The significance of these variables as well as the kinematically allowed inte- 
gration region denoted by q) will be discussed in detail in next section. 

Within the IA the evaluation of (e, e'iVj) cross sections for both proton and neutron knock- 
out determines the inclusive QE cross section. The study of exclusive (e, e'N) reactions 



has been presented in previous work 
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focusing on different aspects of the problem: 
Final-State Interactions (FSI), relativity, correlations, etc. Although such ingredients have 
been proven to be essential in order to fit experimental (e, e'N) cross sections, in what 
follows we restrict our attention to the Plane- Wave Impulse Approximation (PWIA), where 
the knocked-out nucleon has no interaction with the residual nucleus. Being the simplest 
approach to (e, e'N) processes, PWIA retains important relativistic effects that are essential 
in describing reactions at high q and u. Moreover, the (e, e'N) differential cross section in 
PWIA factorizes in two basic terms: the electron-nucleon cross section for a moving, off-shell 
nucleon and the spectral function that gives the combined probability to find a nucleon of 
certain momentum and energy in the nucleus 

da 1 PWIA 
_de'dVt'dp N dVt N \ (e e/Ar) 

with K a kinematical factor [jjj] and where p is the missing momentum and 8 the excitation 
energy, essentially the missing energy minus the separation energy. It is important to point 
out that the factorization property shown in Eq. (121) no longer persists if dynamical relativis- 
tic effects in the bound nucleons are incorporated, i.e., effects from the lower components in 



16Ml8l|. In general we can write 



Ka eN (q, w; p, £, </>n)S(p, £) 



(2) 



the relativistic wave functions, even in the plane- wave limit |20|, |2l| . Note that both the eN 
cross section and the spectral function depend on the two integration variables in Eq. (JTJ), p 
and £. In order to show how the scaling function emerges from PWIA, further assumptions 
are needed. First the spectral function is assumed to be isospin independent, and second 
a eN is assumed to have a very mild dependence on the missing momentum and excitation 
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energy, which is supported by the most commonly used off-shell cross sections Hence 
the eN cross section can be evaluated at fixed values of p and £: typically the differential 
cross section for inclusive QE (e, e') processes is written in the form 

dc 

= a e (q,co;p=\y\,£ = 0)-F(q,u), (3) 

(e,e') 



de'dQ! 



where the single-nucleon cross section is evaluated at the special kinematics p — \y\ (with 
y the scaling variable; see the next section) and £ = (the residual nucleus in its ground 
state). This corresponds to the lowest value of the missing momentum occurring when 
£ = 0. The term a e refers to the azimuthal-angle-averaged single-nucleon cross section and 
it also incorporates the kinematical factor K in Eq. fl2]) and the contribution of all nucleons 



in 



the target, i.e., a e = KJ2t=i I d(f) Ni a eNl /2tt. 



The function F(q, u) in Eq. ((3]) is known as the scaling function and it is given in PWIA 
in terms of the spectral function: 

F{q, w) = 2tt ff pdpd£ S{p, £) . (4) 

A detailed study of the scaling function and its connection with the momentum distri- 
bution will be presented in next section. However, let us start by pointing out some general 
interesting features of this basic result. First, only in the case in which it would be possi- 
ble to extend the kinematically allowed region E(g, oS) to infinity in the excitation energy 
plane, i.e., £ ma x — > oo, would the scaling function be directly linked to the true momentum 
distribution of the A-nuclear system: 

POO 

n(p)= / d£S(p,£). (5) 



J o 

Second, guided by the PWIA result in Eq. ((3]), an experimental scaling function can be also 
defined by dividing the experimental QE (e, e') cross section by the single-nucleon function, 
a e . At high enough values of the momentum transfer q, the function F exp (q,u) has been 

ow the QE peak, that is, F exp becomes only a 



shown to satisfy scaling in the region bei 
function of the scaling variable y (see jl[ 



22| for details). Note that Eq. (0} does 



not apply to F exp (q,u) which incorporates ingredients not included in the simple PWIA 
approach: FSI, MEC, re-scattering processes, etc. The contribution of these effects and 
their impact on the scaling phenomenon depend on the kinematical region explored, leading 
in particular to a significant scaling breaking in the region above the QE peak. 



Furthermore, based on the analysis performed with the Relativistic Fermi Gas (RFG) 
model, and making use of the separate longitudinal (L) and transverse (T) (e, e') data, 
experimental superscaling functions have been introduced: 



fex P (q,u) = k F F exp (q,uj) (6) 
f^P(q,u) = k F F e L x (p( q ,u), (7) 

where hp is the Fermi momentum. In particular, the L response is thought to have very little 
contribution from meson production and from meson-exchange currents and thus should be 
the place where the underlying nuclear dynamics can cleanly be resolved. It has been shown 
to superscale, i.e., the function shows only a very mild dependence upon the momentum 
transfer q (first-kind scaling) and the nuclear system considered (second-kind scaling). This 
has led to introduce a universal experimental superscaling function that constitutes a strong 
constraint for any theoretical model describing QE electron scattering. Not only should the 
superscaling behavior be fulfilled, but also the specific shape of must be reproduced. 
This subject has been studied in detail in previous work showing the importance of FSI and 



relativity |23H27j . and those studies clearly show that any conclusion about the momentum 
distribution based on Eq. (j3J) should be taken with caution. Being aware of this, it is 
illustrative, however, to analyze in detail the basic approaches on which the "link" between 
the momentum distribution and the scaling (superscaling) function is based. Moreover, the 
usual analysis restricted in the past to the region below the QE peak, is now extended to 
the region above the peak, since the superscaling function is defined for both negative 
and positive values of the scaling variable (see discussion in next section). 



II. THE SCALING FUNCTION 



As already shown, in PWIA the scaling function can be expressed as an integral of the 
spectral function S in the (p, £) plane (Eq. (TjJ), with p the struck nucleon's momentum and 

8{p) = yfMg+p* - ^/m° 2 +p* > 0, (8) 

the excitation energy of the recoiling system B, with Mg the ground-state mass of the 
residual nucleus and Mg the general invariant mass of the daughter final state. The inte- 
gration in Eq. (jJJ) is extended to the kinematically allowed region in the (p, £) plane at fixed 
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values of the momentum and energy transfer, (q,ui). This is represented by H(q,oj). The 
genera . 



range 



where 



. kinematics corresponding to QE (e, e') processes leads to the following ^-integration 



max{0,£+} < £ < £~ , 



£ ± (p;q,co) = (M° A + cu) 



M% +p 2 



(9) 



(10) 



(q ±p) 2 + m* N 

and where M\ is the target nuclear mass and the nucleon mass. 

The intercepts between the curve £~ and the p-axis will be denoted by — y and Y, i.e., 
£~(—y;q,uj) = £~(Y;q,cu) = 0. The integration region £(g, oS) is shown in Fig. [T] for fixed 
values of the transferred energy and momentum for u < ujqe (left-hand panel) and cu > ujqe 
(right-hand panel), with uqe the energy where the quasielastic peak (QEP) occurs. In the 
region below the QEP, y is negative and p = —y represents the minimum value for the struck 
nucleon's momentum. Above the QEP y is positive and the curve £ + cuts the integration 
region when p < y. 

q and w fixed ^ q and ro fixed 



(y<o) 




(y>o) 



CO 




FIG. 1: (Color online) Excitation energy corresponding to negative (left) and positive (right) values 
of y. 



In terms of the independent variables q and u the intercepts ±y and Y are given by 



(M° + u)^K 2 - M%W 2 - gA J /W' 
Y(q,u) = + cu)^A 2 - M^W 2 + gA j /W' 



(11) 
(12) 
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with W = \J {M\ + u) 2 — q 2 the center of mass energy and A = (M^ — m 2 N + W 2 )/2. Then 
the scaling function in Eq. (j4j) can be recast as follows 

1 



pdp dSS(p,S) if y< 0(13) 

y Jo 

y r£~(p;q,y) rY(q,y) r£~(p;q,y) 

pdp dSS(p,S)+ pdp I d£S(p,£) if y > (14) 

J£ + (p;q,y) Jy 



'0 J£+{p;q,y) Jy JO 

for negative and positive values of y, respectively. The analysis presented in previous work 
has been restricted to the negative-?/ region, i.e., below the QEP, since this is the region where 
cross section data fulfill y-scaling properties. The function F exp does not scale for positive 
values of y because of the significant scaling violations introduced by effects beyond the IA, 
namely, inelastic processes and contributions from meson-exchange current. However, these 
contributions mostly reside in the purely transverse response and are negligible in the L 
channel. The "universal" superscaling function extracted from the analysis of the separated 
L data, and defined for both negative and positive values of the scaling variable, explains 
our interest in extending the study to the region above the QEP. This strategy, which forces 
us to employ the superscaling function f^ xp to determine F^ xp = f^ xp /k F instead of the usual 
^/-scaling function F exp , can lead to significant effects concerning the momentum and energy 
distribution in the spectral function, as discussed below. 

In the above expressions we have chosen (p, £; q, y) as independent variables. In terms of 
these we can also express the energy transfer 



the limits of the excitation energy 



(q + y) 2 + m 2 N + J M% + y 2 - M' 



A ■ 



yJ(q + y¥ + m? N -yJ(q± P y + m 2 N + M% + y 2 - ^ M% + 



and the upper limit of p: 

Mt (2q + y) + 2(q + y) J ' M% + y 2 y/(q + y) 2 + m 2 N + y [2(q + y) 2 + m 2 N ] 

Y{q,y) = , 

M° 2 + 2y / M° 2 + y ^(q + y y + m 2 N + 2y(q + y) + m 2 N 
In the thermodynamic limit — > oo we get 



(15) 



(16) 



(17) 



£ ± {p- q,y) yj(q + y) 2 + m 2 N - ^ (q ± p) 2 + m 2 N = E q+y - E q±p 
Y(q,y) 2q + y, 



(18) 
(19) 



where we have introduced the nucleon energies = \Jk 2 + m 2 N . Moreover, notice that in 
the limit of very large momentum transfer, i.e., q 3> \y\ and q ^> rrtjv, the above limiting 
values reduce to Y — > 2q and £ ± — > y =F p. 

Following previous arguments presented in , it is instructive to split the spectral 
function into two terms, corresponding to zero and finite excitation energy, respectively: 

S(p,£) = n (p)5(£) + S 1 (p,£) (20) 

with Si(p,£ = 0) = 0, which, inserted in Eqs. (I13II14P yields 

I rY(<i,y) fY(q,y) r£~(p;q,y) 



n(Q,y) rr(q,y) ft {p;q,y> 

—F(q,y<0)= pdpn (p)+ pdp d£S 1 (p,£) (21) 

J_y J_y J 

I rY{q,y) 
—F(q,y>0)= / pdpn (p) 

^ Jv 

i) fY(q,y) f£~{p;q,y) 

S 1 (p,£). 



y r£ (p;q,y) rY(q,y) r£ (p;q,y) 

pdp d£ + pdp d£ 



J£+{p;q,y) Jy JO 



(22) 



In order to analyze how the scaling function and the nucleon momentum distribution are 
connected, we proceed by evaluating the derivatives of the scaling function F with respect 
to y and q. Making use of the Leibniz's formula and choosing (p; q, y) as the three remaining 
independent variables, after some algebra we finally get the following results: 

A. Negative-y region 



if = Y "° (y) O - «"*<-'" + £"* (It) Si(p - r) (23) 
h d i = w(f) + /%*(^)^n. (24) 

Making use of the limits in Eq. (fl6l) and assuming the residual mass to be much 
larger than the momenta, \y\,p, q, we simply have 

d£~ q + y d£~ q + y q — p 



dy E q+y ' dq E q+y E q _ p 
Likewise, the derivatives of Y reduce to dY/dy ~ 1 and dY/dq ~ 2. 



(25) 
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Introducing these results in the general expressions in Eqs. (I23 |24p . we get 



1 dF 
2ir dy 

1 dF 
2^~dq 



Yn (Y) -yn (-y) + 



q + y 



E, 



q+y 



p dp Si{p, £~ 



2Yn Q (Y) 



Y 



pdp 



q+y q-p 



E, 



q+y 



E, 



q-p 



Slip, £~ 



(26) 
(27) 



with £~ and Y given in the thermodynamic limit by Eqs. ( 1181) and (1T91) . Note that the 
excited-state contribution in the spectral function, that is Si, is evaluated at energies along 
the curve £~. 



0.4 



0.3 



> 

O 0.2 



CO 



0.1 



i 1 1 1 j 1 1 1 1 1 1 1 1 1 1 1 r 

'q— >oo 




f 0.2 0.4 0.6 0.8 1 1.2 
-y p [GeV/c] 



1.4 1.6 1.8 



FIG. 2: (Color online) Integration region in the (£,p) plane for y = —0.1 GeV/c and 12 C as the 
target selected. Each curve corresponds to £ ~ for a different momentum transfer. 



For q sufficiently large, q ^> —y, the upper limit Y can be safely taken to oo and, since 
limy^oo Yn (Y) = 0, the expressions for the derivatives simplify to 



q+y J -y 



= ~ ynohy)+ E~ 1 I vdpSxip.S-) (28) 
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FIG. 3: (Color online) Superscaling function f(ip) for negative (a) and positive (b) values of the 
scaling variable ifi. Results correspond to 12 C(e, e') evaluated in RPWIA for different momentum 
transfers. 



If we further assume that Si is small for large values of p, so that the main contribution 
to the integral Eq. (|29|) comes from p ~ — y, then we get 



, dF 
hm — 







(30) 



namely scaling of the first kind (the scaling function F loses its dependence upon q). 

We also observe that, since at a fixed value of y the integration region in Eq. (1271) increases 
with q and the integrand is a positive function, the asymptotic value F(y) is reached from 
below (i.e., monotonically increasing as a function of a) in any PWIA approach, in contrast 



with what experimental data seem to indicate [111, [12|, |22j|. This is clearly illustrated in 
Fig. [2] where the integration region is shown for different values of the momentum transfer 
at fixed y, and it is also consistent with results shown in the left-hand panels (a) of Figs. [3] 
and|H In Fig. [3] we present the superscaling function f(ip) evaluated within the framework 
of the Relativistic Plane- Wave Impulse Approximation (RPWIA) (see [2J, |25( for details) 
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'jO.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.2 0.4 0.6 0.8 



y/k F y/k F 

FIG. 4: (Color online) Superscaling function for negative (a) and positive (b) values of the dimen- 
sionless scaling variable y/kp. Results correspond to 12 C(e,e') evaluated in RPWIA for different 
momentum transfers. 



for different q- values and plotted against the superscaling variable ip in the negative-?/' region 



(below the QEP). This variable is given by 

1 



10 



12| 



X-T 



(31) 



ft*" yj(l + A)r + kv^TTTJ 

where A = u/2m N , n = qj2m^ and r = \Q 2 \/Am 2 N = n 2 — A 2 . The scaling variables y and 
ip are closely connected 121 ] : 



rriir 1 




q 2 2 



y_ 



(32) 



where r/p = kp/m^ and, as noted above, the superscaling function / is connected with F 
via / = kp x F with kp the Fermi momentum. The curves in Fig. E]may be compared with 
the RPWIA results for the superscaling function, now for negative and positive values of 
the dimensionless scaling variable y/kp obtained using the quadratic form of Eq. fl32l ; see 
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Fig. |H As shown, at fixed ip (or y/kp) the function f(ip) increases with q in accordance 
with the previous discussion. The basic results shown in Figs. |3] and H] demonstrate that 
ip and y/kp can be used interchangeably as long as one does not focus on the few percent 
differences seen in the figures, namely, for large magnitudes of the scaling variables. 

In showing the results we choose 12 C as an illustrative example. Indeed this nucleus is 
relevant for many neutrino oscillation experiments, where superscaling ideas can be used to 
make reliable predictions of neutrino- nucleus cross sections 28]. Moreover, the analysis of 
the world data performed in 11] points to an excellent superscaling in the so-called scaling 
region (if; < 0) for nuclei with A > 12. Note, however, that even the 4 He data display a 
very good superscaling behavior for large negative values of the scaling variable (ip < —0.2), 
while at the quasielastic peak there is a 10% violation due to the very different spectral 
function of the lightest nuclei. 



B. Positive-y region 

In this shown in Fig. [1] (right-hand panel), the integration region in the (p,£)- 

plane is limited by the two curves, S + and S~, in the missing momentum region [0, y\. This 
makes the derivative analysis somewhat more complicated. Moreover, the experimental data 
show that scaling arguments of the first kind do not apply to the function F(q,w) in this 
region, that is, F does not become a function only dependent on the scaling variable y. 
On the contrary, it shows a strong dependence upon the momentum transfer q. As already 
mentioned, this is due to important contributions beyond the IA contained in the transverse 
channel. Therefore, although the analysis that follows is applied to F(q,y), it should be 
clearly stated that only the use of the "universal" (namely longitudinal) superscaling function 
/l, and in particular, the study of its derivative with respect to the scaling variable in the 
positive-?/ region, can reveal important effects not accounted for by the results obtained in 
the negative-?/ scaling region. 
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After some algebra, the derivatives of the scaling function F(q, y) are given by 
1 dF fdY\ 

2^ = y "°< y > 



* Y (q,y) f dE~ 

l PdpSfaE-)^ 



- \ pdp Slip, 
'o 



dy 



(33) 



l 9 i= Y ^ Y \ dqJ 



Y{q,y) f dE~ 

pdpSi{p, £~) \ 



y / Qg+ 

pdp Slip, £ + ) f — 



(34) 



As in the previous case, from the general expressions for E ± given in Eq. f[T5"j) and assuming 
the thermodynamic limit, we get 



d£ ± 
dy 

and the derivatives reduce to 



q + y d£ ± q + y q + p 



E, 



q+y 



dq 



E q+y 



E, 



q±p 



(35) 



1 dF 
2tt dy 

1 dF 
2tt 



Yn (Y) - yn (y) + 



q + y 

E, 



q+y 



y 



p dp Si{p, £~ 



p dp Si(p, £~ 



(36) 



2Yn (Y) 



q + y 

E n 



Y 



pdpS^p, £~ 



+ 



q+y Uo 
y q _|_ p 
pdp— — Si{p, £ + 

tiq+y 



pdp Si{p, £~ 



pdp^-^-Siip, £' 

&q+V 



(37) 



Moreover, in the limit of the momentum transfer large enough, q 3> y, so that the 

condition liniy^oo Yu^iY) = holds, the expressions of the derivatives result 

IdF q + y rroo r v 

n — = -y n o{y) + 



2n dy 
1 dF 
2^~dq 



E 



q _|_ y q — p 
pdp [ — ) S^p, £~ 



E 



q+y 



E 



q-p 



y A (q + y q + p\ 

Pdp — — — ) S^p, £ 

V ^q+y ^q+V 



pdpSiip,8 )- / pdpSiip, S + ) (38) 

■)■ 
(39) 

q + y q 

Notice that in the limit in which y can be neglected compared with q, that is, — > —, 

Eg + y Eg 

the same comment applies to the ratio iq+p) / E q+p involved in the second integral in Eq. (13TJD . 
since p is limited within the range [0,y]. Thus, in such a limiting case 

f y pdp (S±l-S±l\siip,£ + )^0 for g»y (40) 

JO V ^q+y ^q+p / 

and only the first integral in Eq. (139]) survives. Furthermore, if the spectral function is such 
that we can neglect p as compared with q inside the integral, we again get scaling of the first 
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kind: lim^oo ^ = 0. This is strictly valid only for very large values of q and it is entirely 
based on the approximations leading to the expression in Eq. (J4]) that connects the scaling 
function with the spectral function. As can be seen from Fig. H] (panel (b), positive-?/ region), 
the RPWIA scaling function shows a negligible dependence on the momentum transfer for 
0-3 < y/kF < 0.8 (q 3> y), whereas for larger y/k F scaling of the first kind begins to be 
slightly violated. The experimental scaling function extracted from the analysis of data at 
intermediate g-values (less than or of the order of the nucleon mass) shows very important 
scaling violations in the region above the QEP (positive values of y). 



1 1 1 ' 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



0.9 
0.8 
0.7 



1 I '/I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 

q=1.0 GeV/c 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 
p [GeV/c] 

FIG. 5: (Color online) As for Fig. [2 but now for positive values of y. Panel (a) corresponds to 
y = 0.1 GeV/c and panel (b) to y = 0.5 GeV/c. 



With regard to the dependence of the scaling function F with q at fixed y, we get different 
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behaviors for small and large values of y. Indeed from Eq. (|39|) we observe that in the case of 
y being very small (in the vicinity of zero), the second integral in Eq. ( 139]) can be neglected. 
As the integrand in the remaining integral is positive, we get dF/dq > 0, i.e., the scaling 
function grows with q. This behavior is in accordance with the one already shown in the 
negative-?/ region. On the contrary, for increasing values of y the first integral in Eq. 
is expected to diminish significantly, since the excitation energy curve £~ along which Si is 
evaluated lies much higher than £ + (see Fig. [5]), and it is reasonable to expect that Si(p, £) 
gets its main contribution for values of the momentum and energy which are not too large. 
For y large enough, only the second integral in Eq. ( 131?j) survives, and as its integrand is also 
positive, the minus sign in front of it leads to dF/dq < 0, that is, the scaling function F 
decreases with q, changing its behavior with respect to the previous cases. It is interesting 
to point out that this result is consistent with the integration regions shown in Fig. [5] where 
for increasing momentum transfer the curve £ + moves to higher excitation energies in the 
{£,p) plane. This means that as q goes up regions at low (£,p) values, where the spectral 
function mostly resides, are not kinematically accessible anymore. A similar argument can 
be applied to the case of very small values of y (see panel (a) in Fig. [5]). However, here the 
integration region lost as £ + goes up with increasing q is less important than the effects 
introduced by the growing integration region attached to £~. This general behavior is also 
in accordance with the RPWIA results for the superscaling function / shown in the right- 
hand panel of Fig. H] (positive values of y), or, alternatively, the right-hand panel of Fig. EJ 
One sees that / increases with q up to ip > 0.4, i.e., yjkp ~ 0.364 (q = 0.5 GeV/c), 
y/kp ~ 0.375 (q = 1.0 GeV/c), y/k F ~ 0.382 (q = oo GeV/c) with k F = 1.2 fm" 1 the Fermi 
momentum. This corresponds to y ~ 0.1 GeV/c, which is the situation represented in panel 
(a) of Fig. |5j Also note that the g-dependence of / in the region where yjkp > 0.3 seen in 
Fig. |4] is very weak. As observed by examining the two panels in Fig. for large y- values the 
energy curves £ ± lie very high, and hence, as q increases, the integrals involved incorporate 
only additional contributions which are very small, leading to a very weak variation with 
momentum transfer. 
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III. NUCLEON MOMENTUM DISTRIBUTION AND THE SCALING FUNC- 
TION 



In the previous section we have derived general integro-differential equations connecting 
the derivatives of the scaling function, dF/dy and dF/dq, with the spectral function. Based 
on these results applied to both negative and positive values of y, in what follows we re- 
visit the "usual" procedure to obtain the nucleon momentum distribution function from the 
analysis of QE (e, e') data. Since the kinematics of electron scattering lead to finite integra- 
tion limits, we may not a priori draw any strong conclusions about the "true" momentum 
distribution given as n{p) = J °° d£ S(p,S), namely the integral of the spectral function up 
to infinite excitation energy. However, assuming the spectral function to reside mostly in 
the (p, £) plane at values of p and £ which are not too large, the previous analyses applied 
to negative- and positive-?/ regions lead to different results, thus providing important and 
complementary information on how the energy and momentum are distributed within the 
spectral function. 



in order to generate the nuclear 



The usual procedure considered in previous work 
momentum distribution from the scaling function has been based on the expression: 

-1 fdFX 



n{k) 



- 



- 



2ny \ dy J 



(41) 



\y\=k 



which has been widely applied in the negative-?/ region. In what follows we extend this 
study to the positive-y region based on the universal superscaling function introduced from 
the analysis of the separated longitudinal data. 

Making use of the general expressions given by Eqs. ( I26|l3"6"j) and assuming the limiting 



case limy^oo Yn (Y) = 0, which is valid if the momentum transfer q is sufficiently large, the 
momentum distribution functions can be written as follows: 



n 



y<0 



n 



y>0 



no(-y) 

n (k) + 

n (y) - 

n (k) - 



_ q + y 

y E q+y 
q-k '°° 



kE q _ k j k 

. q + y 
yEq+y 

q + k 

kE q+k 



pdpSi{p,£~ 
pdpSi{p,£~) 



pdpSi(p,£~ 
pdpS!(p,£~] 



-y=k 



(42) 



p dp Si{p, £ 



pdpSiip,£' 



J y=k 



(43) 
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As observed, both expressions receive contributions from the A — 1 system ground state, 
no(k), as well as from the excited states described through S\(p,£). Although using the 
same notation for the excitation energy £ ~, note that the ^-curves that enter in the spectral 
function Si in Eqs. ( l4"2"j) and (j43p are very different (see Figs. I2"f5"j) . 

Conclusions on the particular behavior of the previous expressions can only be drawn 
based on a specific model for the spectral function; however, it is illustrative to discuss some 
general, "model-independent" , properties. For negative y the function in Eq. (142 p exceeds the 
purely ground-state contribution, i.e., n y<0 (q, k) > no(k) for all q, /c-values. This means that 
the contribution from the excited states adds to the ground-state momentum distribution. 
Concerning the specific role played by each one of the two terms in Eq. (142 p . it is difficult 
to draw stringent conclusions without having control over S±. As the momentum k grows, 
the contribution of the integral in Eq. (142 p is expected to diminish significantly (Si mostly 
residing at momenta and excitation energies which are not too large). A similar comment 
applies also to the ground-state contribution that decreases as k gets larger. The analysis 
of Eq. ( 1431) in the positive-y region differs because of the relative contributions provided by 
the two integrals linked to the excited states. In this case the global response n y>0 (q, k) can 
be smaller and/or larger than the purely ground-state contribution, n (k), depending on the 
specific missing momentum value. 

In what follows we discuss in detail some particular situations, thereby drawing some 
preliminary conclusions on the general behavior shown by n y ^°(q,k). Let us start by con- 
sidering the value of the nucleon momentum k to be in the vicinity of zero. Thus, neglecting 
k as compared with the momentum transfer q (k <C q) and assuming pdp Si(p, £~) 3> 
Jo pdpSi(p,£ + ) — > 0, we can write: 

n y<0 (q,k) ~ n (k) + -%- [ pdp Si(p,£~) > n (k) , (44) 

Jk 

n y >0 (q,k) ~ no (k)-^L J™pdpSi(p,£-)<n (k). (45) 

From these results the following relation (valid for k small enough) occurs, 

n y>0 (q, k) < n (k) < n y<0 (q, k) . (46) 

Moreover, from Eqs. ( f44f45l) the ground-state contribution is roughly given as no(k) ~ 
[n y<0 + n y>0 ]/2. 
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As the nucleoli momentum k grows, the two functions n y<0 (q,k) and n y>0 (q,k) in 
Eqs. ([HI 1431) get closer, crossing each other at some specific k, such that n y>0 (q,k) > 
n y<0 (q, k) for larger k. From the integration region in the (S — p) plane shown in Fig. 0, 
and assuming most of the strength in the spectral function to be located at not too 
high p and S, we can conclude that for intermediate-to-high missing momentum values 
the main contribution in n y>0 (q,k) comes from the second integral in Eq. (1151) . that is, 

n y>°( q ,k)~ 1 £t_f*pdpS 1 (p,£+). 




-1.5 -1 -0.5 0.5 1 1.5 2 2.5 



FIG. 6: (Color online) Average fj^ p (tp) compared with the Gumbel distribution in Eq. (|47h (solid) 
and a fit of the experimental data (dashed). 



To prove these general properties, in what follows we present results based on the deriva- 
tive analysis making use of the superscaling function f(ip). In order to simplify the calcu- 
lations we represent f(ip) by means of the Gumbel probability density function (i.e., the 
derivative of the Gumbel distribution): 



a 



exp 



i$-n) 



a 



exp 



exp 



a 



In our case the values of the parameters are \x = and a = 0.67 (J™ 

cL 



(47) 

/ G (0) = 0.55). 



In Fig. [6] we compare the Gumbel distribution [Eq. ( 1471) ] with fe Xp (ip) and a fit of the 



experimental data 



As shown, the Gumbel distribution nicely fits the data. Moreover, 
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k [fm 1 ] 

FIG. 7: (Color online) Nucleon momentum distribution extracted through the derivative of the 
superscaling function given by the Gumbel probability density in Eq. (|47p . Results corresponding 
to negative (solid line) and positive (dashed) values of the scaling variable are compared. 



it fulfills the unitarity condition f_°° f(ip)dip = 1. The nucleon momentum distribution is 
evaluated through the derivative of the scaling function by using Eq. (14"T]) and recalling that 
/ = k F F, thus getting 



n{k) 



1 1 

27n/ kp dy 

y 



(48) 



J \y\=k 



that, using the approximate relation ip ~ - — , can be presented in the form 



n(k) 



1 1 

2irk kp 



dm 

A k F\^\)\ kFW=k ' 



(49) 



Note that if the superscaling function is not symmetric with respect to ip, as is the case 
for the experimental data, the above expression yields different momentum distributions for 
negative and positive values of ift, which will be denoted by n < and n > , respectively. On the 



contrary symmetric scaling functions, like the RFG one, lead to n < = n > . 
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FIG. 8: (Color online) Integration region in (£ ,p) plane for q = 1 GeV/c and different, negative 
and positive, values of the scaling variable y. The contour curves <? ± in both regions are signaled. 



In the case of the Gumbel distribution we get (setting /x = 0): 

(e-^ - 1) , (50) 



which leads to 



= 2^P ~ X l ^(~*/^) ( 51 ) 



^ = 2^P - e_fc/MF) ] /g( W • (52) 

{ Ayl ( 
In Fig. H we present the results for n^' <0 (k) = — ^ — (solid line) and n^ >0 (k) = — ^ — 

(dashed line) with n^k) and n^{k) given in Eqs. (I51|l52p (at fcj? = 1.2 fm _1 ). As expected, 

rig(fc) and n^{k) (and n^ <0 and n^ >0 , respectively) coincide in the limiting case fc = 0: 

n S (0) = "S (0 ) = ^p- (53 > 

For missing momenta up to k ~ 1 fm" 1 the main contribution resides in n < that is in 
accordance with Eq. (JlSj) and the general discussion presented above. At ~ 1.3 — 1.4 
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fin" 1 , i.e., k close to the Fermi momentum, n < and n > cross each other, with n > being 
much higher for larger k- values. In fact, whereas ?t, < shows a steep slope when k increases, 
which is in accordance with results based on independent-particle model descriptions, n > 
presents a high momentum tail very far from n < and hence from shell-model results (see 
next section). As already explained above, this tail at intermediate-to-high k is linked to the 
much larger contribution given by the spectral function Si when evaluated along the curve 
S + instead of S~. This general behavior is illustrated in Fig. [8] where the contour curves 
£ ± corresponding to positive- and negative-?/ values are presented. The presence of the tail 
at high momentum values in the nucleon momentum distribution is a clear signature of the 
importance of nucleon-nucleon correlations. Since the spectral function maps very different 
regions in the (S — k) plane for negative and positive y (Fig. [S]), the joint analysis of the 
two kinematical regions can provide important clues in the knowledge of NN correlations. 
It should be pointed out that the functions n^ <0 (k) and n^ >0 (k) evaluated through Eq. (149]) 
and presented in Fig. [7] are normalized to different values connected with the different areas 
subtended by the Gumbel distribution function /g(V') a t negative and positive ip, i.e., 0.37 
(for $ < 0) and 0.63 (if> > 0). 



In particular, it has been shown in 
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26] in the framework of relativistic nuclear models 
that the large positive-?/' tail of the scaling function is closely related to final-state interac- 
tions, while the negative-?/' region is more affected by initial-state correlations, as will be 
also shown in the next Section in the CDFM model. The possibility of connecting differ- 
ent aspects of the momentum distribution to initial- and final-state physics will be further 
explored in future work. 

IV. NUCLEON MOMENTUM DISTRIBUTION WITHIN THE COHERENT 
DENSITY FLUCTUATION MODEL 

In this section we give, as an example, the results for the nucleon momentum distribution 



extracted from the scaling function, obtained within the framework o 



model, namely the Coherent Density Fluctuation Model (CDFM) 29 



a particular nuclear 



301 ] . The latter is a 



natural extension to finite nuclei of the relativistic Fermi gas (RFG) model within which 
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the scaling variables ip' was introduced 1 . The CDFM is based on the generator coordinate 



method 



3l| and includes long-range NN correlations (LRC) of collective type. In 



scaling function was defined within the CDFM using the RFG scaling function 



and applied it to various 
In the CDFM model 



processes 



32 



33 



29 



37J4O 



32 
10, 



33 

341 



the 



36] 



30|, the one-body density matrix p(r, r') is an infinite super- 



position of one-body density matrices p x (r,r') corresponding to single Slater determinant 
wave functions of systems of A free nucleons homogeneously distributed in a sphere with 



radius x, density po(x) 



3A 



9tt 



.4 



1/3 



Attx 3 
1.52A 1 / 3 ): 



, and Fermi momentum kp(x) 



3tt 2 



I 1/3 



-po{x) 



a I • 1 
— (with 

X 



00 

p(r,r') = J \F(x)\ 2 p x (r,r')dx. 



(54) 



The weight function |-F(a;)| 2 can be expressed in an equivalent way either by means of the 
density distribution 



29 



30 



33] 



\F(x) 



1 dp{r) 



p (x) dr 



at d -M- < 0, 

dr 



(55) 



or by the nucleon momentum distribution 33] 



\F{x)f 



3tv 2 a dn(k) 



2 x 5 dk 



at < 0. 



k=a/x 



dk 



(56) 



In Eqs. dSSD and (J56 



p(r)dr = A, / n(k)d\t = A, and 



\F(x)\ 2 dx = 1. 



(57) 



In the version of the CDFM approach suggested in 32|, [33], the scaling function has the 
form 



«/(**#'!) 

M)= f \F(x)\ 2 f RFG (x,ij')dx, 
'0 



(58) 



1 The scaling variable ijj' differs from ip by a phenomenological energy shift E s ~ 20 MeV (for 12 C) intro- 
duced in order to reproduce the experimental position of the quasielastic peak: ip'(q,uj) = i^ili w — E s ). 
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where the RFG scaling function is 



kpx\ip' 



a 



xm,N\ 2 (kpx\i\)' 



a 



a 



x 



a 



xm N 



-2 A 1 + 



xm N 



(59) 



In the CDFM the Fermi momentum kp is calculated for each nucleus by 



k, 



k F (x)\F(x)\ 2 dx = / -\F(x)\ 2 dx 



a , 



x 



(60) 



and is not a fitting parameter, as it is in the RFG model. 

By using Eqs. ( 1551) and ( |56i) in Eqs. ( |58i) and ( 160|) . the CDFM scaling function / (■*/>') and 



fcp can be expressed equivalently by the density and momentum distributions 

a/{kwW\) 



33] 



4:71 

~A 



p{x) 



x 2 f RFG (ip' , x) + 



where f RFG (ip',x) is given by Eq. (159|) . and 



4tt 



n{k) 



k 2 f RFG ^',k) + 



x 3 df RFG (ip' , x) 
3 dx 



k 3 df RFG ^',k) 
3 dk 



dx, 



(61) 



(62) 



k F W\ 



where 



/rfg(^',&) = ^ 



m N y ( k F \ip'\ f ^ 



x 



2 + 



2W1 + 



m N 



• (63) 



Eq. ( 1621) is valid under the condition 



lim [n(k)k 3 ] = 0. 



(64) 



k— >oo 

^From Eq. ( 162]) one can estimate the possibility to obtain information about the nucleon 
momentum distribution from the empirical data for the scaling function. If we keep only 
the main term of the RFG scaling function from Eq. ( 1631 : 

2" 



/rfg ty, k) 



1 



kptjj' 
~~k~ 



(65) 
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and its derivative 



then 



2 



dk 



- I n(k)k 2 

'k F \ip'\ 



k 3 



3 k 2 



dk. 



In Eq. Q£ 



n(k)dk = 1. 

Using Eq. (EZJ), n{k) can be found by solving the integral-differential equation: 

oo 

i dfW) 



(66) 



(67) 



(68) 



n(k) 



2nk 2 d(k F \ij'\) 



k F \tp'\=k 



dk'n(k') 



(69) 



In this work we solve the above equation from CDFM using the experimentally obtained 
scaling function. The latter can be represented the Gumbel probability density function in 
Eq. (I47p . The results for the nucleon momentum distribution obtained in this way are given 
in Fig. [9] by dashed lines in both cases: n < (/c) for ip < (green dashed line) and n > (fc) for 
i/j > (blue dashed line) . They are compared with the results obtained using the expression 
for n(k) through the derivative of the scaling function, Eq. (j49|) . 

The momentum distributions n < {k) and n > {k) obtained by using Eq. (I49p and the ex- 
perimental scaling function presented by Eq. (1471) are given in Fig. [9] by solid lines. For a 
comparison we present in the same figure the momentum distributions from the RFG model 
(7Zrfg) ; the shell-model results (using Woods-Saxon single-particle wave functions) for 56 Fe 
(%s)i as we h as the momentum distribution (t^lfd) obtained within the Light-Front Dy- 
namics (LFD) approach 4l| (see also 33] and the late modification of the approach in 38]). 
The latter is based on the nucleon momentum distribution in the deuteron (including its 
high-momentum component) from the LFD method {e.g., 42|,|43[ and references therein). In 
the calculations kp = 1.2 fm _1 . In Fig. [9] all nucleon momentum distributions are normalized 
to unity [Eq. (|68l)]. 

One can see from Fig. [9] that, in general, the results for n(k) in CDFM confirm the 
considerations made in Sects. IhMlHI Namely, 

(i) At k < 1.3 fm -1 the CDFM momentum distributions (from Eq. ( 1691 ) ra < (/c) > n > (/c), 
while at k > 1.3 fm _1 n < (/c) < n > (/c). The same is valid for the momentum distribution 
obtained using Eq. P9|) . This is in accord with the general consideration from Sect. IHIl 
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FIG. 9: (Color online) The nucleon momentum distribution extracted from the scaling function. 
Solid lines: n < (light lines, green online) and n > (dark lines, blue online) obtained through the 
derivative of the scaling function [Eq. (|49p ]: dashed lines: n < (light lines, green online) and n > 
(dark lines, blue online) using the CDFM integral-differential equation [Eq. (I69j) ]. The Gumbel 
probability density function /g(V ; ) [Eq. (I47p ] is used in the calculations. For comparison are given 
the momentum distributions from the Relativistic Fermi Gas model (wrfg)) from the shell model 
( n ws) an d from the Light Front Dynamics (tilfd)- All momentum distributions are normalized to 
unity [Eq. (pi) ]. 

(ii) The crossing point of the lines showing n < (/c) and n > (/c) from Eq. (I49p is at slightly 
smaller value of k than that for n < (/c) and n > (/c) obtained from Eq. (!69|) . This follows 
also from the comparison of the explicit forms of Eqs. fl69|) and (l49|) . 

(iii) n < (/c) from Eq. ( |69i) is close to the result for n < (fc) from Eq. ( I49l . while the difference 
between Tv'ik) from Eq. (H9|) and n > (/c) from Eq. fl69|) increases with k. At k = 4 fm -1 
rv'ik) from Eq. (149|) is around twice larger than n > (/c) from Eq. (169|) . At the same 
time for fc < 1.2 fm" 1 , n > (/c) from Eq. fl69l) is larger than Tv'ik) from Eq. (jUJ). 
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V. CONCLUSIONS 



In the present work a study of the scaling function and its connection with the mo- 
mentum distribution is presented. As is well known, a close relationship between the two 
quantities exists using the PWIA and under some conditions for the kinematically allowed 
region (E(q,oj)), once one has accounted for the roles of FSI, MEC, rescattering processes, 
etc. Here these restricted approximations are considered in detail. The "usual" analyses 
performed in the past to the region below the QE peak is extended to the region above 
the peak, since the superscaling function is defined for both negative and positive values 
of the scaling variable. This is justified, since a "universal" superscaling function has been 
extracted from the analysis of the separated longitudinal data. The explicit expressions for 
the derivatives dF/dy and dF/dq for both negative- and positive-;?/ regions are derived and 
their dependences on q and y are analyzed. 

The general integro-differential equations connecting the derivatives dF/dy and dF/dq 
with the spectral function are derived. The results obtained allow us to revisit the "usual" 
procedure to obtain the nucleon momentum distribution from the analyses of the QE scat- 
tering data. The considerations in the present work lead to results that are quite different 
from those obtained solely in the negative-?/ scaling region and give information about the 
energy and momentum distribution in the spectral function. It is shown that the expressions 
for the nucleon momentum distributions n y<0 (q,k) and n y>0 (q,k) have contributions from 
the momentum distribution no(k) of the ground state of the system with A — 1 nucleons, 
as well as from the part of the spectral function Si(p,£) that contains information about 
the excited states. It is shown that for small momenta k: n y>0 (q,k) ^ n Q (k) ^ n y<0 (q,k), 
while as k grows the two functions n y<0 (q,k) and n y>0 (q,k) get closer, crossing each other 
at some value of k and yielding n y>0 (q, k) > n y<0 (q, k) for higher k. 

The general properties of the momentum distribution established in the present work are 
validated by the results obtained from the derivative analysis using the superscaling function 
f(ip) represented by the parameterized Gumbel probability density function that provides a 
good fit to the experimental longitudinal scaling function f^p{ip)- It is concluded that the 
high-momentum tail of the momentum distribution is a clear signature for the important 
effects stemming from nucleon-nucleon correlations. 

The general properties of the nucleon momentum distribution obtained are also illustrated 



26 



using the scaling function obtained in the framework of a particular nuclear model, namely 
the Coherent Density Fluctuation Model (CDFM) that includes collective long-range NN 
correlations. It is shown that the momentum distribution in the CDFM has the properties 
already pointed out in the general consideration. 
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